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ABSTRACT 

We investigate the non-linear evolution of the matter power spectrum by using a large 
set of high-resolution N-body/hydrodynamic simulations. The linear matter power in 
the initial conditions is consistently modified to mimic the presence of warm dark 
matter particles which induce a small scale cut-off in the power as compared to stan- 
dard cold dark matter scenarios. The impact of such thermal relics is examined at 
small scales k > 1 /iMpc~^, at redshifts of z < 5, which are particularly important for 
the next generation of Lyman-a forest, weak lensing and galaxy clustering surveys. 
We measure the mass and redshift dependence of the warm dark matter non-linear 
matter power and provide a fitting formula which is accurate at the ^ 2% level below 
z = 3 and for particle masses of hiwdm^ 0.5 keV. The role of baryonic physics on the 
warm dark matter induced suppression is also quantified. In particular, we examine 
the effects of cooling, star formation and feedback from strong galactic winds. Finally, 
we find that a modified version of the halo model describes the shape of the warm 
dark matter suppressed power spectra better than HALOFIT. In the case of weak lensing 

however, the latter works better than the former, since it is more accurate on the relevant, mid-range 
scales, albeit very inaccurate on the smallest scales (fc > lO/iMpc"^) of the matter power spectrum. 

Key words: Cosmology: theory - large-scale structure of the Universe - dark matter, 
methods: numerical - gravitational lensing: weak 
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1 INTRODUCTION 

The increasing amount of observational data available and 
the numerical tools developed for their interpretation have 
given rise the so-called era of precision cosmology. At the 
present time, the cosmological concordance model based on 
a mixture of cold dark matter and a cosmological constant 
must thereby be tested in new regimes (both in space and 
in time) and using as many observations and techniques as 
possible in order to either confirm or disprove it. 

Among the many different observables the non-linear 
matter power spectrum is a crucial ingredient since it allows 
us to describe the clustering properties of matter at small 
scales and low-redshift, where linear theory is not reliable. 
However, accurate modelling of non-linear physical processes 
is needed, in order to use this observable to gain quantitative 
results on the nature of dark matter. 

Warm Dark Matter (WDM) is an intriguing possibil- 
ity for a dark matter candidate. Its velocity dispersions 
are intermediate between those of cold dark matter and 



hot dark matter (e.g. light neutrinos). In this scenario, at 
scales smaller than WDM free-streaming scales, cosmolog- 
ical perturbations are erased and gravitational clustering 
is significantly suppressed. If such particles are initially in 
thermal equilibrium, they have a smaller temperature and 
affect smaller scales than those affected by neutrinos. In ad- 
dition, WDM produces a distinctive suppression feature at 
such scales as compared to that induced by neutrinos. For 
example, thermal relics of masses at around 1 keV which 
constitute all of the dark matter have a free-streaming scale 
that is comparable to that of galaxies, well into the non- 
linear regime. Among the different WDM candidates a spe- 
cial role is played by the s t erile ne utrino with mass at the 
keV scale (|Bovarskv et al.l (|2009al )). WDM was originally 
proposed to solve some putative problems that are present 
in cold dark matter scenar ios at small scales (see lColin et al.l 
(|200Cl ): lBode et al.l ||200H )). however it is at present contro- 
versial whether these tensions with cold dark matter pre- 
dictions can be solved by modifying the nature of dark 
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matter particles or by som e other baryonic process (e.g. 
iTruiillo-Gomez et all |20ld )). 

In the present paper we wish to quantify the im- 
pact of a WDM relic on the non-linear power spec- 
trum by using a set of N-body/hydrodynamic simulations 
of cosmological volumes at high resolution. Investigating 
WDM scenarios in a cosmological setting has been done 
by means of N-body codes in order to carefully quan- 
tify the impact of such a candidate in terms of halo 
mass function, structu r e formation, halo d ensity proper- 
ties (iBode et al.l (|200ll ): IColm et al.l (|2008l ): IColombi et aH 
lJ2009f )^ and particular care needs to be pla ced on correctly 
addre ssing numerical/convergence issues (jWang fc White! 
l|2007f )l. In general, while the WDM induced suppression 
transfer func tion can b e relia bl y estimated in th e linea r 
regime (e.g. I Viel et al.' ("20051); iBovanovskv et al.1 (|2008l ): 
iLesgourgues fc Tram (2011)), the non- linear suppression has 
not been investigated. However, a recent attempt to obtain 
the non-linear matter power at small scales by modify ing 
the halo model is described in lSmith fc Markovid (|201ll '). 

The analysis of matter power spectra at small scales 
has been performed in recent year by different groups 
by focussing mos tly on bary o n phy s ics such as f e edbac k 
and cooling ( e.g. 'Rudd ct al.' ('2008'); 'Guillet et all (|2010l ): 
ICasarini et al . (2011); van Daalcn ct al. (2011)). 

Recently, Ivan Daalen et al.l (J2011i ) presented an ex- 
tensive investigation of the effects of several different 
implementations of galactic feedback on the total mat- 
ter power. Including feedback from Active Galactic Nu- 
clei (AGN) was claimed to be most realistic since it 
matches optical and X-ray observations of groups of 
galaxies and s olves the overcooling problem (see also re - 
lated works bylPuchwein et al.1 (12003'): iFabianet al.l (|201ol ): 
iMcCarthv et al.l (|2010l '): iTevssier et al.l (|201ll )')~This sce- 
nario has a relatively large impact in terms of total mat- 
ter power at the scales affected by the presence of a WDM 
candidate. In a subsequent paper, ISemboloni et al.l (|201in 
carefully analysed the impact of this feedback on weak lens- 
ing observables and concluded that it will not be possible 
to constrain WDM masses with future weak lensing surveys 
(as claimed in for example iMarkovic et al.l (|201ll )V These 
findings are very interesting and show that a future mea- 
surement of the matter power at such small scales should 
be considered with an accurate model of baryonic physics 
and not only of the dark matter component. In the present 
work, although we will be exploring some feedback mecha- 
nism that do not include the AGN feedback, the main focus 
is on the signature of the WDM suppression in terms of to- 
tal matter power, presented as a ratio with a corresponding 
ACDM standard model that includes the same astrophysical 
input and differs only in the initial total matter power. 

Different constraints can be obtained by using sev- 
eral astrophysical probes. For example by using Lyman- 
Q observables such as the transmitted Lyman- a flux 
power, very competitive measurements in the form of 
lower limits (mwDM > 4 keV, 2a C.L.) have been de- 
rived by using the SDSS flux powe r and other high e r red- 
shift and higher resolution data (|Viel et al.1 |2005| . l2006l : 
ISeliak et al.l 120061 ): these constraints become much weaker 
if the WDM is assumed to account to only a given frac- 
tion of the dark matter (JBovarskv et al.l l2009al ) or if the 
initial linear suppression for a sterile neutrino is consid- 



ered (jBovarskv et al.l l2009bl ) (in this latter case basically 
any mwDM, sterile > 1 keV is allowed). Alternatively, con- 
straints on WDM models can be placed using: the evo- 
lution and size of small scale stru cture in the local vo l- 
ume high resolution simulations (jTikhonov et al.l 120091 '): 
simulated M ilky Way haloes to probe properties of sate l- 
lite galaxies (|Polisenskv fc RicottilbOllHLovell et al.llJOllI '): 
large scale structure data (| AbazaiianI 12006) : the formation 
of th e first stars and gal axies in high resolution simula- 
tions (jGao fc Theung|2007); weak lensing power spectra an d 
cross-spectra ( Markovic et al.l 201ll: Semboloni et aLlbOllh : 
the dynamics of the satel lites (Knebe et al.ll2008l V the abun- 
dance of sub-structures (Co lin et al. 200d): the inner prop- 



erties of dwarf galaxies (Strigari et al.ll2006l '): mass function 
in t he local group as d etermined from radio observations in 
HI (jZavala et alj|2009|): the clus tering properties of galaxies 
at small scales ( Coil et al.ll2008l ): the properties of satellites 
as inferred from semi-anal ytical models of galaxy formation 
(|Macci6 fc FontanotlpOlQ): phase-space densit y constraints 
from dwarf galaxies (|de 'Vega fc Sancheaboid ). 

We believe that most of the astrophysical probes used 
so far in order to constrain the small scale properties 
of dark matter could benefit from a comprehensive nu- 
merical modelling of the non-linear matter power. The 
present work aims at providing such a quantity by using N- 
body/hydrodynamic simulations. The findings could also be 
useful for future surveys such as PanSTARRS, DBS, LSST, 
ADEPT, EUCLID, JDEM or eROSITA, WFXT and SPT. 

The layout of the paper is as follows. In Section 2 we 
present our set of simulations and the code we use in order 
to investigate the non-linear suppression on the total mat- 
ter power. Section 3 contains the main results of the present 
work and the description of the checks made in order to 
present a reliable estimate of the WDM non-linear suppres- 
sion: we focus on numerical convergence, box-size, baryonic 
physics, particle velocities and the effect induced by cosmo- 
logical parameters on the WDM power. As an application of 
the findings in Section 3 we present the weak lensing power 
and cross-spectra for a realistic future weak lensing survey 
in Section 4 and compare these results with those that could 
be obtained by using either linear-theory or halo models in 
Section 5. We conclude with a summary in Section 6. 



2 THE SIMULATIONS 

Our set of simulations has been run with the parallel hydro- 
dynamic (TreeSPH: Tree-Smoothed Particle Hydrodynam- 
ics) code GADGET-2 based on the conservative 'entropy- 
formulation' of SPH (|Springej[2005l '). Most of the runs use 
the TreePM (Tree-Particle Mesh) N-body set-up and con- 
sist only of dark matter particles, however for a few runs, 
in order to test the impact of baryonic physics, we switched 
hydrodynamic processes on. 

The cosmological reference model corresponds to a 
'fiducial' ACDM Universe with the following parameters, at 
z = 0, ^m = 0.2711, f^A = 0.7289, fib = 0.0451, n^ = 0.966, 
Ho = 70.3 km s"^ Mpc"^ and ag. = 0.809. This model is in 
agreement w ith the recent constra ints obtained by WMAP- 
7 year data (JKomatsu et al.ll2011h and by other large scale 
structure probes. The initial (linear) power spectrum is gen- 
erated at 2 = 99 with the publicly available software CAMB 
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Table 1. Summary of the simulations performed. Linear box-size, 
mass of warm dark matter particle and gravitational softening 
are reported in comoving units (left, center and right columns, 
respectively). The particle- mesh (PM) grid is chosen to be equal 
to Nj-,j^ with NjjM = 512^. Simulations (a) have been run with 
hydrodynamic processes (a simplified star formation recipe and 
radiative processes for the gas) and with full hydrodynamics with 
the standard multiphase modelling of the interstellar medium and 
strong kinetic feedback in the form of galactic winds. Simulations 
(a) have been also run at lower resolution A^dm = 384"^ and for 
different values of o-g , Om and Hi). Simulation (b) has been run 
by switching the initial velocities of warm dark matter particles 
off and by increasing the linear size of the PM grid by a factor 
3. Simulation (c) has been run with Nqm = 640"^ dark matter 
particles with a softening of 1 kpc/h to z = 0.5. 



U and then modified to simulate warm dark matter (see be- 
low). 

We consider different box sizes in order to address both 
the large scale power and (more importantly) the effect of 
resolution. The gravitational softening is set to be 1/40-th of 
the mean linear inter-particle separation and is kept fixed in 
comoving units. The dimension of the PM grid, which is used 
for the long-range force computation, is chosen to be equal 
to the number of particles in all but a single case in which a 
finer grid is used. The simulations follow a cosmological peri- 
odic volume filled with 512^ dark matter particles (an equal 
number of gas particles is used for the hydrodynamic simu- 
lations) in all but two cases in which a smaller and a larger 
number of particles are chosen in order to check for numeri- 
cal convergence of matter power. We mainly focus on WDM 
masses around 1 keV. For such a mass, the characteristic cut- 
off in the power spectrum appears at scales of about k '^ 1.5 
ftMpc"^, the suppression reaching 50% at fe = 6 /iMpc~^. 
The suppressed scales are highly non-linear and therefore re- 
quire high-resolution as well as N-body techniques. However, 
in order to be conservative we present results for the fol- 
lowing mwDM values: 0.25, 0.5, 1, 2 and 4 keV. These limits 
could be easily converted to masses for a sterile neutrino par- 
ticle produced in the so-called standard Dodelson-Widrow 
scenario and correspond to rus = 0.7, 1.66, 4.4, 11.1, 28.1 keV 
(note that physically motivated scenarios based on for ex- 
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ample non-resonant production mechanisms have been pro- 
posed, however the simulations carried out in the present 
work cannot be strictly applied to those since they require 
a non-trivial modification of the line ar transfer function, as 
discussed bv lBovarskv et all (|2009bl )). 

The initial conditions for warm dark matte r particles 
are ge nerated using the procedure described in IViel et al.l 
(2005), which we briefly summarize here. The linear ACDM 
power is multiplied by the following function: 

'5/v 



TiL(fc) 



PwDM(fc)/PACDM(fc) = (1 + {aky 

IkeV \^" /^v— -°" 



Q(mwDM) = 0.049 ( 



mwDM 



/f2wDM\ ■ f ^ \ ( 

\ 0.25 J Vary ^ 



where v = 1.12 and a has units of h"^ Mpc (e.g. 
[Hansen et al.l 120021 ) . We stress that the above equation is 
an approximation which is strictly valid only at fc < 5 — 10 
h Mpc~^ . Below this scale the warm dark matter power spec- 
trum could be described by a more complicated function and 
acoustic oscillations are present (se e for example the recent 
work in iLesgourgues fc Tramlboill ). 

Initial velocities for warm dark matter particles are 
drawn from a Fermi-Dirac distribution and added to the 
proper velocity assigned by linear theory: the r.m.s. veloc- 
ity dispersion associated to their thermal motion is 27.9, 
11.5, 4.4. 1.7, 0.7 km/s for mwDM=0.25,0.5,l,2,4 keV, re- 
spectively. The typical r.m.s. velocity dispersion for the dark 
matter particles of the ACDM runs is ~ 27 km/s, so at least 
for masses above 1 keV the thermal WDM motion is a small 
fraction of the physical velocity dispersion assigned by the 
Zel'dovich approximation. 

When baryonic physics is included, we consider the fol- 
lowing processes: i) radiative cooling and heating, ii) star 
formation processes. Hi) feedback by galactic winds. We 
note that metal cooling is not included in the simulations 
and only cooling from H and He is considered, while galac- 
tic winds are powered by massive stars and not AGN. 

The rationale is to see at which level these processes 
impact the non-linear matter power at small scales in order 
to compare to the suppression effects of WDM. Thus, we are 
not aiming at exploring in a comprehensive way the impact 
of t hese processes on the non-linear power at sm all scales, 
(e.g. Ivan Daalen et al.ll201ll : ICasarini et al.ll201ll ): the bary- 
onic simulations are used only to quantify the impact of such 
processes on the suppression induced by WDM w.r.t. cold 
dark matter scenarios. 

Radiative cooling (H and He) as well as heating pro- 
cesses are assumed for a primordial mix of hydrogen and 
helium corresponding to a mean Ultraviolet Background 
similar to th at produced by q uasars and galaxies and im- 
plemented in lKatz et al.l IjlQQo ). This background naturally 
gives a hydrogen ionization rate r_i2 ~ 1 at high redshift 
and an evolution of the physical state of the intergalactic 
medium (IGM) wh ich is in agreement with observations (e.g. 
iBolton et al.ll2005h . The star formation criterion for the de- 
fault runs is a very simple one that converts all the gas 
particles whose temperature falls below 10^ K and whose 
density contrast is larger than 1000 into coUisionless stars 
(more details can be found in IViel et al.l 12004 ). This pre- 
scription is usually called "QLYA" (quick Lyman-a ) since 
it is very efficient in quantitatively describing the Lyman- 
Q forest and the low density IGM. We also run a simula- 
tion with the full multi-phase description of the interstel- 



4 M. Viel et al. 



ACDM 



WDM 




-1.0 -0.5 0.0 0.5 1.0 1.5 2.0 

log (h-5d„) 

Figure 1. "Visual" inspection of the redsiiift evolution of cosmic structures in the ACDM and WDM (mwDM=l keV) scenarios (left 
and right columns, respectively) for the defaults (25,512) runs. From the top to the bottom rows we show a 2.5 h~^ Mpc thick slice of 
the projected dark matter density at z = 0, 2, 5 respectively. At z = the clustering properties of the dark matter at scales fc < 10 
/iMpc~-'^are indistinguishable in the two scenarios, while at z = 2,5 the WDM model has a suppression in power of about 5% and 25% 
at fc = 10 hMpc-l. 
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lar medium (ISM) and with kinetic feedback in the form 
of strong galactic winds as in ISpringel fc Hernguistl (12003) ■ 
The chosen speed of the wind is 483 km/s and both the 
ISM modeUing and this feedback mechanism is expected to 
impact on the distribution of baryons and thus on the to- 
tal matter power spectrum. The wind particles temporarily 
decouple from the hydrodynamics: the maximum allowed 
time of the decoupling is fdec = I/ v-k with 1=20 kpc/fe and 
vw_= 483 km/s (see for e xamp le iDalla Vecchia &: Schavd 
l|2008h and lPierleoni et all (|2008l ') to understand the effects 
of this parameter in terms of feedback efficiency and on the 
properties of local Hl-galaxies). 

We note that simulations that include baryons are sig- 
nificantly slower than the default dark matter only runs and 
therefore our constraints will mainly be derived from the 
former simulations. 

In the following, the different simulations will be de- 
noted by two numbers, (A^i, N2): N\ is the size of the box in 
comoving Mpc//i and N2 is the cubic root of the total num- 
ber of gas particles in the simulation. The mass per dark 
matter particle is 8.7 x lO^'Mo/ft for the default (25,512) 
simulations. This mass resolution allows to adequately sam- 
ple the free-streaming mass for the models considered here. 

In Figure [T] we show the projected dark matter density 
as extracted from the default (25,512) runs in the ACDM 
case (left) and WDM case (right) for mwDM = l keV. This 
WDM particle mass is already ruled out at a significant level 
by Ly man- Q forest observations (e.g. ISeliakll2005l : IViel et al.l 
120061 ). The different rows refer to 2: = 0, 2, 5 from top to 
bottom, respectively. In this Figure it is essential to see how 
the clustering proceeds differently in the two scenarios and 
while there are large differences below the Mpc scale at z = 5 
between the two cosmic webs, these differences are largely 
erased by non-linear evolution at 2 = 0, 2. 

The main features of the simulations are summarized in 
Table 1. 



3 RESULTS 

In this section we describe the main results obtained from 
our sample of simulations. The power is computed from the 
distributions of the different sets of particles (dark matter, 
gas and stars) separately and for the total matter compo- 
nent by performing a CIC (Cloud-In-Cell) assignment to 
a grid of the size of the PM grid. The CIC kernel is also 
deco nvolved when g etting the density at the grid points 
(e.g. IViel et al.ll20ld )). We also show a small scale estimate 
(fe > 10 /iMpc~^) of t he power obtained with the fold- 
ing m ethod described in (jJenkins et al.ll 19981 : IColombi et al.l 
I2OO91 ). although this power will not be used quantitatively. 
We will plot the suppression in power as a percentage 
difference between WDM and ACDM total matter power 
spectra, normalized by the default ACDM total matter 
power. The initial conditions for CDM and WDM have the 
same phases and cosmological/astrophysical parameters in 
order to highlight the effect of the warm dark matter free 
streaming. 



3.1 Resolution and box-size 

In Figure [2] we show the percentage difference between 
the total non-linear powers of WDM (iuwdm^I keV) and 
ACDM runs. We subtract the shot-noise power from all the 
power spectrum estimates made. For our largest box-sizes 
the shot-noise power is comparable to the actual measured 
power at z = at fc ~ 150 /iMpc~^, while for the de- 
fault simulations (25,512) of iuwdm^I (0.25) keV the matter 
power is always above the shot-noise level for 2; < 10 and for 
fc<20(7) /iMpc-\ 

This figure focusses on the resolution and box-size ef- 
fects and presents the percentage difference at four different 
redshifts z = 0, 1, 3, 5 (bottom right, bottom left, top right 
and top left panels respectively) and for three different box- 
sizes (100, 50, 25 h~^ Mpc shown as blue, black and green 
curves respectively). The dotted line represents the redshift 
independent linear cut-off of Eq[l] while the lower resolu- 
tion (25,384) run is also plotted in orange. Here, there are 
two estimates for the power: one at large scales (continu- 
ous curves), the second at smaller scales (dashed curves). 
We are primarily interested in the power at scales fc < 10 
/iMpc~^and therefore only the large scale estimate will be 
used. However, we also show the power at smaller scales 
since physical and numerical effects play a larger role in 
this range. We note that the linear theory suppression is 
a good approximation only at fc < 1 /iMpc^^. From the 
figure one can see that there is convergence up to fc = 50 
/iMpc"^between (25,512) and (25,384) runs in the redshift 
range considered. The resolution used is thus suffici e nt for 
mwDM=l keV particles. Note that Ivan Daalen et al.l (|2011h 
recently found that (100,512) ACDM simulations have suf- 
ficiently converged at scales fc < 10 /iMpc~^. At fc = 3(10) 
h~^ Mpc and z — ^ there is already a 5 (50)% difference 
between the linear and non-linear power. At 2: = 0, 1, 3 the 
differences between WDM and ACDM power is below 1%, 
2% and 5% respectively at fc = 10 /iMpc~^. The maxi- 
mum suppression dip is strongly influenced by resolution 
and moves to larger wavenumbers when the resolution in- 
creases. At fc > 100 /iMpc~^we note a steep (resolution de- 
pendent) turn-over in the suppression which is likely to be 
due to effects that impact on the halo structure and which 
has also been found in CDM numerical simulations that in- 
clude a fr action of the matter content in the form of active 
neutrinos (JBrandbvge et al.ll2008l : Iviel et al.ll2010h . 

We have checked that increasing the particle-mesh grid 
by a factor three (i.e. PM=1536) has negligible impact on 
the total matter power at scales fc < 100 /iMpc~^. In order 
to test the robustness of our results in terms of shot-noise 
level we have also run a WDM simulation with iuwdm^I 
keV and A'^dm = 640'^ particles and compared the power 
spectra with the (25,512) and (25,384) runs: we confirm 
very good agreement between these simulations at fc < 20 
h Mpc~^in the redshift range considered in the present work. 
More precisely, the (25,512) and (25,640) WDM runs agree 
below the one percent level at fc < 100 /iMpc^^. 



3.2 The effect of the mass of a warm dark matter 
particle 

Here we address the effect of varying uiwdm on the non- 
linear matter power. The results are shown in Figure [3] 
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Figure 2. Percentage difference between warm dark matter non-linear power and cold dark matter for the different runs. Tlie mass of 
the warm dark matter particle is kept fixed to mwDM = 1 keV. Blue, black, green curves refer to 100, 50, 25 h~^ Mpc respectively and 
with a fixed number of particles JVdm = 512'^. The orange curves refer to 25 h~^ Mpc and has a fixed number of particles A^dm = 384'^ 
The continuous lines represent the large scale estimate of the power, while the dashed ones describe the small scale power obtained with 
the folding method (see text). The four panels represent different redshifts at z = 0, 1, 2, 5 (bottom right, bottom left, top right and top 
left, respectively). The dotted line plotted at z = and 2 = 5 is the redshift independent linear suppression between the two models. 



where we report five different masses for tlie (25,512) de- 
fault runs. Tlie curves correspond to niwDM=0.25,0.5,l,2 
and 4 keV (orange, green, black, blue and red curves, re- 
spectively) at 2; = 0, 1,2 and 5 (bottom right, bottom left, 
top right and top left, respectively). The linear suppressions 
are also shown with dotted lines of the corresponding colors. 
At 2 = 5 we can see large differences between the models 
that become smaller with the redshift evolution. The 20% 
suppression at fc = 10 /iMpc~^at 2 = 5 for the mwDM = l 
keV model becomes 2% at 2 = 1 and it is below 1% at 2 = 0: 
the clustering properties of the dark matter are the same at 
scales above fc ~ 10 /iMpc~^at least for mwDM> 1 keV. 
The mwDM=0-5 keV model still presents a 7% suppression 
by 2 = 0, while the suppression is four times larger at 2 = 2. 
The linear suppression is a very poor approximation in the 
range of wavenumbers considered here even at high redshift. 
At 2 = 1, which is particularly interesting for weak lensing 



data, a 2% measurement of the non-linear power is likely to 
be able to exclude models below the 1 keV value (bottom 
left panel). The dip of the maximum suppression and the 
turn-over both move to larger scales as the mass decreases. 
We have also investigated the importance of WDM ve- 
locities in the initial conditions by running a simulation 
without assigning a Fermi-Dirac drawn thermal velocity to 
the dark matter particles. We tested this for a mwDM = l 
keV model and found differences always below 1% in terms 
of total matter power at the scales of interest here. 



3.3 Baryonic effects 

In this section we explore the effects of baryonic physics 
on the warm dark matter suppression. Baryons amount to 
about 17% of the total matter content and we expect that 
astrophysical processes affecting their properties can impact 
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Figure 3. Percentage difference between warm dark matter non-linear power and cold dark matter for the different runs. The resolution 
is kept fixed in this plot and only 25 h~^ Mpc boxes are considered. Orange, green, black, blue and red curves refer to mwDM = 
0.25, 0.5, 1, 2, 4 keV, respectively. The continuous lines represent the large scale estimate of the power, while the dashed ones describe the 
small scale power obtained with the folding method (see text). The four panels represent different redshifts at ^ = 0, 1, 2, 5 (bottom right, 
bottom left, top right and top left, respectively). The dotted coloured curves plotted at z = and z = 5 are the redshift independent 
linear suppression between the different models. 



the total matter power at small scales at some level. We 
identify three processes that are able to modify the clus- 
tering properties of baryons: radiative processes, star for- 
mation and galactic feedback. These processes are usually 
modelled by hydrodynamic simulations of galaxy formation. 
Here, the main goal is not to explore fully the many param- 
eters governing these important physical aspects, but rather 
to address their impact in WDM models by adopting pre- 
scriptions that are widely used in the literature. There could 
well be other astrophysical processes (radiative transfer ef- 
fects, feedback from active galactic nuclei, etc.) that can 
also affect the distribution of baryons and their c lustering 
properties (see for example lvan Daalen et al.ll201ll ). 

In Figure [4] we plot the WDM suppression for the de- 
fault simulation of mwDM=l keV for three different cases: 
pure dark matter (green curves); a hydrodynamic simula- 
tion that includes cooling (H and He only) and heating 



by an ultraviolet background as well as the simple star 
formation criterion able to simulate the Lyman-a forest 
("BARYONS+QLYA" run in blue); a hydrodynamic sim- 
ulation that includes the full star formation model based on 
the multi-phase description of the ISM (sophisticated com- 
pared to QLYA) and strong galactic feedback in the form 
of winds. ("BARYONS+SF+WINDS" in black). Unfortu- 
nately, due to the fact that hydrodynamic simulations are 
slower than dark matter only runs it was not possible to 
carry this last simulation down to 2 = and it was stopped 
at 2 = 1.2. 



We also report the ratio between mass in stars 
f^WDM/"lcDM(^ = 0,1,2,3) = 0.85,0.78,0.7,0.6 
for the "BARYONS-HQLYA" runs, while we have 
^wdm/^acdmI-z: = 1-2,3) = 0.87,0.63 for the 
"BARYONS-fSF+WINDS"runs with galactic winds feed- 
back. For both the WDM and ACDM runs the mass fraction 



© 0000 RAS, MNRAS 000, 000-000 



8 M. Viel et al. 



I 

3 
< 

a. 



« 
« 



-10 
-20 
-30 
-40 
-50 




DM only 

BARYONS+QLYA 
BARYONS+SF+WINDS ^ 



10 



3 



« 
« 











1 


1 ' / - 






^^\\ 


/ ' : 
/ ; : 


-10 


r 




% 


; ; : 
\ ; ' - 


-20 


L 




x; 


\y-> ; ; 




i Z- 


=3 




\\ ; 


-30 


'r 






U ;' \ 


-40 


r 






\ \ / '-- 


-SO 









). 1 



100 



10 




10 




^^'^^^"^ 




: 


20 






\ 1 

\ ''^ 


,! 


30 


; Z-- 


=1.2 


\ \ \ 

\\ \ 

1 


\ , - 


40 


L 




\ 
\ 


; - 








\ 


^ / : 


50 




1 




: 



100 




1 10 100 

k (h/Mpc) 



1 10 100 

k (h/Mpc) 



Figure 4. Percentage difference between total WDM non-linear matter power and total ACDM non-linear matter power for runs that 
incorporate baryonic physical processes. The simulations refer to a 25 h~^ Mpc box and mwDM=l keV. The green curves refer to the 
pure dark matter simulations; blue curves refer to simulations that include baryons, cooling from H and He and a simplified recipe for star 
formation that turns into coUisionless stars all the gas particle below T=10 K and denser than 1000 times the mean density (QLYA); 
black curves are instead obtained by using the default criterion of multi-phase star formation of Springel (2005) and feedback in the form of 
strong kinetic driven winds (this simulation was stopped at 2 = 1.2). The continuous lines represent the large scale estimate of the power, 
while the dashed ones describe the small scale power obtained with the folding method (see text). The four panels represent different 
redshifts at 2 = 0, 1.2, 3, 5 (bottom right, bottom left, top right and top left, respectively). In the z = panel (note the different scale for 
the J/— axis) we also show as the red and cyan curves the percentage difference between total matter power spectra that include and do 
not include cooling for ACDM (red) and WDM (cyan) models; namely the quantity 100 X (pbary°ns+QLYA _ pDMONLYj/pDMONLY 



in stars is reduced by a factor five when winds are included 
compared to the "BARYONS+QLYA" case. 

All of these processes can significantly change the clus- 
tering of baryons especially at intermediate scales where 
baryon pressure is important (fc ~ 1 /iMpc~^), where they 
are not expected to trace the dark matter and at smaller 
scales due to the complex interplay between feedback and 
star formation processes. Cooling as well as heating modify 
the thermal properties of the gas and are important espe- 
cially for the low density IGM; the star formation criterion 
determines how much gas is turned into stars within the po- 
tential wells of dark matter haloes; galactic winds displace 
gas out of the galaxies into the low density IGM, usually 



in a hot phase that prevents subsequent cooling. Since the 
cosmic structure is generally different in CDM and WDM 
models we do not expect the WDM suppression to be ex- 
actly the same between two simulations that share the same 
astrophysical prescriptions. From Figure |3J one can see that 
dark matter only simulations are in good agreement (at the 
percent level up to fc = 10 hMpc~^) with simulations that 
include radiative cooling (metal cooling is not included) and 
QLYA star formation, while at smaller scales there are signif- 
icant differences. It is clear that the presence of baryons and 
star formation greatly affects the maximum suppression and 
the turn-over. Note that differences much larger than 10% 
between simulations implementing different radiative pro- 
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cesses (e.g. metal cooling) or feedback recipes a,re expected 
at k > 20 fcMpc'Hn ACDM models (see e.g. iRudd et all 
I2OO8I : iGuillet et aLllJOlol : Ivan Daalen et al.ll201lh . Further^ 
more, in the case of AGN feedback at the level required to 
match observed gas fractions of groups a 10% difference is 
found already at k = 1 /tMpc~^and a 1% re duction already 
at fc = 0.3 feMpc~^ (|van Daalen et all (|201ll ')'). 

In the z = Q panel we also show the difference in the 
power spectra of ACDM and WDM models by normaliz- 
ing to the corresponding dark matter only model, in or- 
der to highlight the effect of cooling produced by baryons 
as opposed to the WDM signature. The two percentage 
differences are shown as cyan (WDM) and red (ACDM) 
curves: the WDM universe when filled with baryons that 
can cool has more power than a corresponding ACDM 
universe filled with the same baryon fraction. The quan- 
tity Pni,wDM,cooiing/-Pni,WDM,dmoniy is about 5% larger than 

^nl, ACDM, cooling/-Pnl, ACDM, dmonly at fc = 10 /l MpC~ and 

z = 5, at z = 1 it becomes only 2 % larger and by 2: = 0, 
there are no differences between the two quantities at A: = 10 
ftMpc"^. The cooling of baryons inside the potential wells 
of dark matter haloes produces further collapse of structures 
and in general increases the (total) matter power spectrum. 
It is thus likely than in the WDM model the baryons cool 
slightly more efficiently than in the corresponding ACDM 
since at high redshifts, the collapse of haloes around the 
WDM cutoff is rapid and small scale modes affected by cool- 
ing (H and He) gro w more rapidly t h an in CDM: this is also 
the trend found by iGao fc Theund (|2007l ) from the analy- 
sis of cooling at very high resolution and high redshift in 
hydrodynamic simulations. 

The WDM suppression is thereby highly infiuenced by 
astrophysical effects at A: = 100 /iMpc^^. In general we ex- 
pect an additional suppression due to baryons of about 2-3% 
at A: = 10 /iMpc~^at z > 1.5 for mwDM = l keV, while this 
discrepancy becomes smaller at lower redshifts. The num- 
bers above do not apply once AGN feedback is included and 
are greatly underestimated, if AG N feedback impa c ts the 
matter power at the level found bv lvan Daalen et al.l ()201ll ) 
and Scmboloni ct al. (2011). In order to accurately measure 
power on these scales, any such AGN feedback should be 
accounted for. 



3.4 Other cosmological parameters 

To test the robustness of our results we extended the set 
of simulations by exploring also other cosmological parame- 
ters, namely: Q,-m, Ho and erg. In order to do that we modify 
the input linear ACDM parameter calculated by CAMB and 
vary one parameter at a time. It is clear that some param- 
eters like (Tg (or As) do not have any impact at the linear 
level, while they could impact the non-linear power in a 
way that should by quantified with simulations. We choose 
the following parameters for the WDM and correspond- 
ing ACDM runs: fim = 0.22, 0.32, Ho = 62, 78 km/s/Mpc 
and (Tg = 0.75, 0.87. When calculating the suppression we 
always normalize both simulations to the same ag value 
(erg = 0.75,0.809 and 0.87). Since the WDM suppression 
of the power spectrum has a relatively distinct shape and 
the cut-off scale is at much higher k than what is probed by 
erg normalisation, we expect WDM to be nearly independent 
of any other parameter probed. The range explored by the 



Ho values produces a maximum ±2% difference in terms of 
the WDM suppression compared to the reference Ho = 70.3 
km/s/Mpc case ai k — 1 — 10 h Mpc~'^and at z < 3, while at 
2: = 5 there is a 5% difference at A; = 10 /iMpc"^. The fini 
parameter produces a maximum difference of 1% at z < 3 
in the same range of wavenumbers and about 5% at 2 = 5 
and A; = 10 /iMpc~^. A different choice of erg has a slightly 
larger impact. This is seen in Figure [5]), where the WDM 
induced suppression with such a choice of erg is divided by 
the reference case of erg = 0.809. It is clear from the figure 
that the large (10 %) differences in place at z = 5 are largely 
canceled by the non-linear growth and are at the ±2% level 
at 2 = 1 — 2 and at the 3% level at A = 10 ^Mpc~^today. 

Motivated by the present findings we regard our non- 
linear cutoff and its redshift dependence as robust at least for 
the range of cosmological parameters investigated at z < 3, 
for iuwdmSS 0.5 keV and at A; = 1 — 10 /iMpc^^: in fact 
the differences are at the ± 2% level and in the next section 
we will provide a fitting formula with a comparable level of 
accuracy. Larger masses for uiwdm will only result in smaller 
differences in terms of WDM suppression. 

We also notice that degenerate features with the non- 
linear WDM suppression might arise in the context of non- 
standard models of da rk energy, a s e.g. interacting dark en- 
ergy scenarios (see e.g. lBaldilboid ). The investigation of such 
possible degeneracies goes beyond the scope of the present 
paper. 



3.5 An analytical fitting formula 

Inspired by the corresponding formula for the linear sup- 
pression, we have found the following fitting formula to be 
a good approximation of the late time evolution of the non- 
linear suppression with an accuracy at the 2% level at z < 3 
and for masses larger than mwDM=0.5 keV: 

Tii{k) = PwDM{k)/PACDM{k) = (1 + {aky')-"^", 

IkeV \^-^^ n + z-^ " 



t(mwDM, 2) = 0.0476 ( 



(2) 



- mvVDM / 

with 1/ = 3, I = 0.6 and s — 0.4. 

We have chosen as a pivot redshift 2: = 1 since this is the 
redshift where accurate weak lensing data will be available. 
This formula has been derived from the dark matter only 
runs. 



WEAK LENSING SHEAR POWER 
SPECTRA 



Following iMarkovic et al.l (|201ll ) and ISmith fc Markovid 
([2011j) , we examine the effect of the fitting function in Equa- 
tion [2] on the weak lensing power spectrum. Weak gravita- 
tional lensing is the distortion found in images of distant 
galaxies due to the deflection of light from these gala:x- 
ies by the gravitational potential wells of intervening mat- 
ter. F or a review, see for example iBartelmann fc Schneideij 
(J200l|). The advantage of gravitational lensing is that un- 
like other large scale structure data, it does not require a 
knowledge of galaxy bias for the derivation of the prop- 
erties of the underlying dark matter density field and is, 
at least on large scales, independent of baryonic physics. 
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Figure 5. Impact of a different erg value in terms of WDM-induced suppression. The four panels represent different redshifts at z = 
0, 1, 2, 5 (bottom right, bottom left, top right and top left, respectively) for the (25,512) with mwDM=l keV. Green represents the 
{erg = 0.809) reference case, while the two other curves indicate the suppression for ag = 0.87 (black) and eg = 0.75 (blue). 



In other words, the weak lensing power spectrum directly 
probes the matter power spectrum. However, weak lensing 
measures the matter power spectrum at low redshifts. For 
this reason it is necessary to have available robust models 
of non-linear structure. For a survey able to probe angular 
multipoles from Z ~ 20 up to i ~ 2 x 10*, in the redshift 
range of 2 = 0.5 — 2.0, the corresponding range of wavenum- 
bers must be fc ~ 0.005— 15 /iMpc^^. Note that the matter 
power at fc > 10 /iMpc~^only has a significant contribution 
to the weak lensing power spectrum at lower redshifts, where 
however the lensing power is lower. 

Future weak lensing surveys accompanied by extensive 
photometric redshift surveys will be able to disentangle the 
contribution to weak lensing by dark matter at different red- 
shifts , by binning source galaxies into tomographic bins (jHul 
1 1993 ). By cross and auto correlating the lensing power in 
these bins, the three dimensional dark matter distribution 
can be reconstructed. An existin g example of such a recon- 
struction is the COSMOS field (|Massev et al.ll2007l ). Such 



tomography probes the non-linear matter power spectrum 

at different redshifts. 

We use HALOFIT (jSmith et al.ll2003l) to calculate non- 
linear corr ections to the approximate linear matter power 
spectrum (|Malll996l ). We then apply Equation [5] to approx- 
imate the W DM effects and find t he weak lensing power 
spectrum (e.g. IXakada fc Jainll2004l ): 



a,{i) 



dxm{Xi)Wjixi)XiPni ( fc = -,Xi 



Xi 



(3) 



where Xi(^i) is the comoving distance to the lens at redshift 
zi and Wi is the lensing weight in the tomographic bin i: 



w/ ^ 47rG' 

ai(zi)c^ 



ni[zB) J — r—dzB , (4) 



where we assume a fiat universe and ai{zi) is the scale factor 
at the redshift of the lens, pm.o is the matter energy density 
today and ni{zs) is the normalised redshift distribution of 
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Figure 6. The percentage WDM effect in auto- and cross-correlation power spectra of redshift bins at approximately 2 = 1 and z = 1.6, 
respectively. All the lines are calculated from non-linear matter power spectra modified for WDM by the fitting function in Equation 
[2] for WDM particle masses of 1 keV (left panel) and 0.5 keV (right panel). In addition we plot predicted error bars for a future weak 
lensing survey, dividing the multipoles into 20 redshift bins. Note that the error bars on auto and cross power spectra of different bins 
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black line at 0) models, one must know the entire covariance matrix for a survey. Note secondly that the auto power spectra of redshift 
bins at 2 = 1 and z = 1.6 have an upturn around I ~ 10^. This is due to the dominance of shot noise on those scales. This upturn is not 
present in the cross power spectrum, because through cross correlation this noise due to intrinsic galaxy ellipticities is eliminated. 



sources in the i-th tomographic bin. We bin the multipoles 
into 20 bins. 

In order to assess detectability of WDM by future weak 
lensing surveys, we calculate predicted error bars on the 
weak lensi ng power spectrum us ing the covariance matrix 
formalism (jTakada fc JainI 120041 ') and assu ming errors for 
a futu re re alistic weak lensing surve y as in iMarkovic et al.l 
(|201l[) and lSmith fc Markovi? (120111 ) with 8 redshift bins in 
the range z = 0.5 — 2.0. We plot the resulting percentage 
differences between WDM and CDM weak lensing power 
spectra in Figure [SI It is important to note that the error 
bars in the figure do not fully characterise the sensitivity of 
the power spectra, since there are additional correlations be- 
tween the error bars of different bin combinations. Addition- 
ally, there are correlations in the error bars on large I (small 
scales) due to non-linearities. Further statistical tests using 
the entire covariance matrix must be used in order to fully 
account for the above correlations. For this plot we choose 
only the 5— f/i and 8— i/i redshift bins, whose source galaxy 
distributions have the mean at z ~ 1.0 and 1.6 respectively. 
These bins are chosen because they represent a range with 
the maximal WDM effect as well as lensing signal. Note that 
the upturn around I ~ 10^ in the auto-correlation power 
spectra of bins 5 and 8 is due to the dominance of shot 
noise on those scales. This noise is due to intrinsic galaxy 
ellipticities and can be eliminated by cross-correlating dif- 
ferent redshift bins, as can also be seen in Figure [6] (see also 
iTakadafc j"airill2004l ). 

In the right panel of Figure [S] we plot the effects of 
the 0.5 keV particle and since the black dashed line lies 
far outside the error bars this is a strong indication that 



such a particle can be ruled out (or detected) by a future 
weak lensing survey. This is consistent with pre vious works 
(|Markovic et al.ll201ll : ISmith fc Markovidl201ll ). In the left 
panel of Figure [6] we plot the effects of a 1 keV WDM parti- 
cle: in this case it is more difficult to distinguish from CDM 
(black dashed line), but the strongly affected cross power 
spectra are still significantly different from their expected 
values in ACDM. 

In a recent paper. ISemboloni et al.l (|201ll ) have explored 
the effect that AGN feedback has in terms of matter power 
and weak lensing power spectra finding that there is a sup- 
pression of about 30% at A: = 10 ^Mpc~^when this feed- 
back mechanism is included. Although they did not investi- 
gate WDM models, this result is important, since it shows 
that the effect could be much larger than the correspond- 
ing WDM induced suppression and comparable at 2 = 
to the mwDM=0.25 keV case. It is clear that future weak 
lensing surveys aiming at measuring the matter power at 
these scales should carefully consider AGN effects since they 
could be degenerate with cosmological parameters such as 
the mass of the WDM particle. 



5 COMPARISON WITH HALO MODEL 

As described in Section [H it is necessary to have a robust 
model of non-linear structure in order to take full advan- 
tage of future weak lensing data. For this reason we com- 
pare the non-linear matter power spectra extracted from 
our simulations with previously derived non-linear models. 
The halo model of non-linear structure is based on the 
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Figure 7. The comparison of different non-linear models at redshifts 1.0 (top panels) and 0.5 (bottom panels) for WDM particles with 
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simulations from previous plots with the fiducial values for trg- The blue solid lines are the corresponding analytical fits from equation [2] 
The green solid lines are calculated using the modified halo model, whereas the green dashed line is the standard HALOFIT. The dotted 
line is the effect as seen in the linear matter power spectrum. 



assumption that large scale structure is made up of indi- 
vidu al objects occupying peaks in the matter overdensity 
field (|Press fc Schechtedll974l : ISeliakll2000l : ICooray fc ShethI 
l2002h . The most important elements of this model , the mass 
function, the halo b ias (IPre ss & Schccht_eiJ [l97J) and the 
halo density profile (jNavarro et al.lll995l ) are based on the 
assumptions that all dark matter in the universe is found in 
haloes and that there is no observable suppression of small 
scale overdensities from early-times free-streaming of dark 
matter particles or late-times thermal velocities. 

These are characteristic pro perties of CDM, but do not 
apply to WDM. For this reason ISmith fc Markovid (120111 ) 
modified the halo model by applying a specific prescription 
to the non-linear contribution, in addition to suppressing the 
initi al density field, m odelled by applying a transfer function 
from I Viel et al.l ((20051) to the linear matter power spectrum. 
Such prescription consists of: i) treating the dark matter 
density field as made up of two components: a smooth, linear 
component and a non- linear component, both with power at 
all scales; ii) introducing a cut-off mass scale, below which 
no haloes are found; Hi) suppressing the mass function also 
above the cut-off scale and iv) suppressing the centers of 
halo density profiles by convolving them with a Gaussian 



function, whose width depended on the WDM relic thermal 
velocity. 

Here, we do not attempt to explore each of these el- 
ements with simulations individually, but rather compare 
the final matter power spectra found fro m simulations with 
those from the WDM halo model of ISmith fc Markovig 
(l201ll) . 

Secondly, ISmith et al.l (|2003l ) compared the standard 
CDM halo model to CDM simulations of large scale struc- 
ture formation and developed an analytical fit to the non- 
linear corrections of the matter power spectrum, known as 
HALOFIT. We apply the se corrections to a linear matter 
power suppressed by the IViel et al.l (|2005l ) WDM transfer 
function (see Equation [T]) . 

We show the results of these comparisons in Figure 
[7] As before, we plot the percent differences between the 
WDM and CDM matter power spectra obtained from our 
simulations of WDM only. We show this for particle masses 
of mwDM=l keV (left panels) and mwDM=0.5 keV (right 
panels) at redshifts z = 1 (top row) and z = 0.5 (bot- 
tom row). We find that the WDM halo model is closest 
to simulations at redshift 1 for 1 keV WDM, but that it 
over-estimates the suppression effect at redshift 0.5 for 0.5 
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from cross power spectra of the 5th and 8th tomographic bins (corresponding to z = 1.6 and z = 1, respectively) for WDM particle 
masses of mwDM=l keV (left panel) and mwDM=0.5 keV (right panel). 



keV WDM by about 5 percent on scales fc > 1. On scales 
fc < 1 /iMpc~^however, the HALOFIT non-linear correction 
describes the simulations better than the halo model, even 
though on smaller scales it severely underestimates the sup- 
pression effect, which becomes worse at lower redshifts. A 
further small modification of the WDM halo model will im- 
prove its correspondence to the simulations and allow one 
to use it at small scales. 



We additionally consider these models of non-linear 
WDM structure to calculate the weak lensing power spec- 
tra in order to explore the significance of using the correct 
model. We again plot percentage differences between WDM 
and CDM weak lensing power spectra in Figure [8l We show 
only curves representing the cross correlation power spec- 
trum of redshift bins at z = 1 and z — 1.6 for consistency 
with Figure [6l We again examine WDM models with parti- 
cle masses of mwDM = l keV (left panel) and mwDM=0.5 keV 
(right panel). We also calculate the weak lensing power spec- 
tra without non-linear corrections to the matter power spec- 
trum and note that this severely over-estimates the effect of 
WDM suppression. In the lensing calculation, the HALOFIT 
non-linear corrections applied to the WDM suppressed lin- 
ear matter power spectrum seem to perform better in de- 
scribing the results of our WDM simulations than than the 
WDM halo model. This due to the fact that the range of 
wavenumbers that are better described by the HALOFIT cor- 
rections, namely k < lhMpc~^ are significantly more rele- 
vant to the weak lensing power spectrum than the smaller 
scales where HALOFIT strongly deviates from the simulation 
results. 



6 CONCLUSIONS 

By using a large set of N-body and hydrodynamic simula- 
tions we have explored the non-linear evolution of the total 
matter power. The focus of the present work is on small 
scales and relatively low redshifts where non-linear effects 
are important and need to be properly modelled with sim- 
ulations. We checked for numerical convergence and box- 
sizes/resolution effects in the range fc = 1 — 10 hMpc~^ . We 
explored how different masses of a warm dark matter candi- 
date affect the non-linear suppression as compared to a cor- 
responding ACDM model that shares the same parameters 
and astrophysical inputs. Our findings can be summarized 
as follows: 

- Cosmological volumes of linear size 25h''^ comoving 
Mpc and with 512"^ DM particles are sufficient to sample 
the WDM suppression for nivvDM^ 1 keV at the percent 
level at fc < 10 /iMpc"\ 

- The non-linear suppression induced by WDM is 
strongly redshift dependent. However, by z = 0, up to fc = 10 
hMpc~^ , there are virtually no differences (below 1%) be- 
tween ACDM and WDM models with mwDMSS 1 keV. 

- At higher redshifts differences are larger, being closer 
to the linear suppression. At 2; ~ 1 there are differences of 
the order of a few percent between the non-linear WDM and 
ACDM power spectra. 

- Baryonic physics and in particular radiative processes in 
the gas component, the star formation criterion and galactic 
feedback in the form of winds are likely to affect the matter 
power at the 2-3 % level in the range fc = 1 — 10 hMpc~^. 
However, a much stronger e ffect can be expe c ted if AGN 
feedback is considered as in (|van Daalen et al.l (|201ll )). 
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- We investigate how a change in Sim, ^o and erg im- 
pacts the non-linear power and WDM suppression in par- 
ticular, when values different from our reference choice are 
used. Small difference are found (at the ± 2% level) at the 
scales considered here. Thus the WDM cutoff has a distinc- 
tive feature which is not degenerate with other cosmological 
parameters also at a non-linear level. 

- We provide a useful fit to the non-linear WDM induced 
suppression in terms of a redshift-dependent transfer func- 
tion; this fitting formula should agree to the actual measured 
power at the 2% level at z < 3 and for masses above 0.5 keV. 

- Reaching a higher accuracy (percent level) in terms of 
WDM non-linear power would require a very extensive anal- 
ysis of astrophysical aspects related to the baryonic compo- 
nent such as considering different feedback effects. Among 
these the most promising seems to be AGN feedback which 
happen to solve the overcooling problem and has a strong 



impac t on the total m a tter p ower (see Ivan Daalen et al 
(|201ll ): ISemboloniet all l|201ll )'). 

- We find that future weak lensing surveys will most 
likely be powerful enough to rule out WDM masses smaller 
than 1 keV, which i s con sistent with previo u s res ults of 
iMarkovic etaP (J201ll l and ISmith fc Markovij (J201ll ). Rul- 
ing out models for masses larger than 1 keV would still be 
possible by using the cross-correlation signal between differ- 
ent redshift bins. However, measurement of the weak lensing 
power at these scales should also consider the effect of bary- 
onic physics car efully and parameters could be biased as 
recently found in lSemboloni et al.l l|201ll ). where it has been 
shown that AGN feedback produces a suppression which is 
larger than the one induced by WDM at the scales consid- 
ered here. 

- Non-linear corrections to the matter power spectrum 
in the WDM scenario obtained from HALOFIT correspond 
better to the results of the WDM only simulation at scales 
A; < 10 /iMpc~^, if compared to the non-linear corrections 
of the halo model from Smith & Markovic (2011). Because 
these scales are most relevant for weak lensing power spec- 
tra, using HALOFIT yields a better correspondence to the 
weak lensing power spectra calculated using our fitting func- 
tion. However, on scales fc > 10 /iMpc~^, the halo model 
performs slightly better in that it better describes the shape 
of the suppression in the power spectrum, even if it does 
overestimate the effect. For this reason we believe that a 
further modification to the halo model may be needed, es- 
pecially for weak lensing power spectra calculations. 



As recently show n by Ivan Daalen et al.l (|201l|) and 
ISemboloni et al.l (J201ll ) including AGN feedback has strong 
consequences in terms of matter power and weak lensing, 
a comprehensive analysis that aims at measuring the mass 
of a warm dark matter candidate should thus hope to lift 
the degeneracies present (i.e. suppression in terms of matter 
power) by exploiting the different redshift and scale depen- 
dencies and by fully exploring the astrophysical parameter 
space and marginalize over the nuisance parameters. 

We believe that future efforts aiming at measuring the 
coldness of cold dark matter should investigate the non- 
linear matter power in the range z = — 5 either using weak 
lensing observables or the small scale clustering of galaxies. 
These constraints can be particularly useful since they are 
complementary to those that can be obtained from high red- 



shift Lyman-a forest data (e.g. BOSS/SDSS-IH survey) or 
galactic and sub-galactic observables in the local universe. 



ACKNOWLEDGMENTS. 

We are grateful to Robert Smith for providing tables to 
calculate the halo model power spectrum and to the ref- 
eree Joop Schaye for his comments. Numerical computations 
were performed at the High Performance Computer Cluster 
Darwin (HPCF) in Cambridge (UK) and at the RZG com- 
puting center in Garching. We acknowledge support from 
an ASI/AAE grant, AST contracts Euclid-IC (1/031/10/0), 
INFN PD51, PRIN MIUR, PRIN INAF 2009, ERC-StG 
"cosmoIGM", the DFG Cluster of Excellence "Origin and 
Structure of the Universe", the TRR33 Transregio Col- 
laborative Research Network on the "Dark Universe", and 
from the European Commissions FP7 Marie Curie Initial 
Training Network CosmoComp (PITN-GA-2009-238356). 
Post-processing of the simulations has also been carried 
at CINECA (Italy) and COSMOS Supercomputer in Cam- 
bridge (UK). KM acknowledges support from the Interna- 
tional Majc-Planck Research School. 



REFERENCES 

Abazajian K., 2006, Phys. Rev. D, 73, 063513 

Baldi M., 2010, MNRAS, in press 

Bartelmann M., Schneider P., 2001, Physics Reports, 340, 

291 
Bode P., Ostriker J. P., Turok N., 2001, ApJ, 556, 93 
Bolton J. S., Haehneh M. C, Viel M., Springel V., 2005, 

MNRAS, 357, 1178 
Boyanovsky D., de Vega H. J., Sanchez N. G., 2008, Phys. 

Rev. D, 78, 063546 
Boyarsky A., Lesgourgues J., Ruchayskiy O., Viel M., 

2009a, JCAP, 5, 12 
Boyarsky A., Lesgourgues J., Ruchayskiy O., Viel M., 

2009b, Physical Review Letters, 102, 201304 
Brandbyge J., Hannestad S., Haugb0lle T., Thomsen B., 

2008, JCAP, 8, 20 
Casarini L., La Vacca G., Amendola L., Bonometto S. A., 

Maccio A. V., 2011, JCAP, 3, 26 
Coil A. L., Newman J. A., Croton D., Cooper M. C, Davis 

M., Faber S. M., Gerke B. F., Koo D. C, Padmanabhan 

N., Wechsler R. H., Weiner B. J., 2008, ApJ, 672, 153 
Colin P., Avila-Reese V., Valenzuela O., 2000, ApJ, 542, 

622 
Coh'n P., Valenzuela O., Avila-Reese V., 2008, ApJ, 673, 

203 
Colombi S., Jaffe A., Novikov D., Pichon C, 2009, MN- 
RAS, 393, 511 
Cooray A., Sheth R., 2002, Physics Reports, 372, 1 
Dalla Vecchia C, Schaye J., 2008, MNRAS, 387, 1431 
de Vega H. J., Sanchez N. G., 2010, MNRAS, 404, 885 
Fabjan D., Borgani S., Tornatore L., Saro A., Murante G., 

Dolag K., 2010, MNRAS, 401, 1670 
Gao L., Theuns T., 2007, Science, 317, 1527 
Guillet T., Teyssier R., Colombi S., 2010, MNRAS, 405, 

525 



© 0000 RAS, MNRAS 000, 000-000 



Non-linear Matter Power in WDM Cosmologies 15 



Hansen S. H., Lesgourgues J., Pastor S., Silk J., 2002, MN- 

RAS, 333, 544 
Hu W., 1999, The Astrophysical Journal, 522, L21 
Jenkins A., Frenk C. S., Pearce F. R., Thomas P. A., Col- 
berg J. M., White S. D. M., Couchman H. M. P., Peacock 
J. A., Efstathiou G., Nelson A. H., 1998, ApJ, 499, 20 
Katz N., Weinberg D. H., Hernquist L., 1996, ApJS, 105, 

19 
Knebe A., Arnold B., Power C, Gibson B. K., 2008, MN- 

RAS, 386, 1029 
Komatsu E., et al., 2011, ApJS, 192, 18 
Lesgourgues J., Tram T., 2011, ArXiv e-prints 
LoveU M., Eke V., Frenk C., Gao L., Jenkins A., Theuns 
T., Wang J., Boyarsky A., Ruchayskiy O., 2011, ArXiv 
e-prints 
Ma C.-R, 1996, Astrophysical Journal v.471, 471, 13 
Maccio A. V., Fontanot F., 2010, MNRAS, 404, L16 
Markovic K., Bridle S., Slosar A., Weller J., 2011, JCAP, 

1, 22 
Massey R., Rhodes J., Leauthaud A., Capak P., Ellis R., 
Koekemoer A., Refregier A., Scoville N., Taylor J. E., Al- 
bert J., Berge J., Heymans C., Johnston D., Kneib J. -P., 
Mellier Y., Mobasher B., Semboloni E., Shopbell P., Tasca 
L., Waerbeke L. V., 2007, arXiv, astro-ph 
McCarthy I. G., Schaye J., Ponman T. J., Bower R. G., 
Booth C. M., DaUa Vecchia C, Grain R. A., Springel V., 
Theuns T., Wiersma R. P. C, 2010, MNRAS, 406, 822 
Navarro J. F., Fronk C. S., White S. D. M., 1995, arXiv, 

astro-ph 
Pierleoni M. 

282 
Polisensky E., Ricotti M., 2011, Phys. Rev. D, 83, 
Press W. H., Schechter P., 1974, ApJ, 187, 425 
Puchwein E., Sijacki D., Springel V., 2008, ApJ, 687, L53 
Rudd D. H., Zentner A. R., Kravtsov A. V., 2008, ApJ, 

672, 19 
Seljak U., 2000, Monthly Notices of the Royal Astronomical 

Society, 318, 203 
Seljak U., Makarov A., McDonald P., Trac H., 2006, Phys- 
ical Review Letters, 97, 191303 
Seljak U. e. a., 2005, Phys. Rev. D, 71, 103515 
Semboloni E., Hoekstra H., Schaye J., van Daalen M. P., 

McCarthy I. G., 2011, ArXiv e-prints 
Smith R. E., Markovic K., 2011, ArXiv e-prints 
Smith R. E., Peacock J. A., Jenkins A., White S. D. M., 
Frenk C. S., Pearce F. R., Thomas P. A., Efstathiou G., 
Couchman H. M. P., 2003, MNRAS, 341, 1311 
Springel V., 2005, MNRAS, 364, 1105 
Springel V., Hernquist L., 2003, MNRAS, 339, 289 
Strigari L. E., Bullock J. S., Kaplinghat M., Kravtsov 
A. v., Gnedin O. Y., Abazajian K., Klypin A. A., 2006, 
ApJ, 652, 306 
Takada M., Jain B., 2004, Monthly Notices of the Royal 

Astronomical Society, 348, 897 
Teyssier R., Moore B., Martizzi D., Dubois Y., Mayer L., 

2011, MNRAS, 414, 195 
Tikhonov A. v., Gottlober S., Yepes G., Hoffman Y., 2009, 

MNRAS, 399, 1611 
Trujillo- Gomez S., Klypin A., Primack J., Romanowsky 

A. J., 2010, ArXiv e-prints 
van Daalen M. P., Schaye J., Booth C. M., Dalla Vecchia 
C, 2011, ArXiv e-prints 



Branchini E., Viel M., 2008, MNRAS, 388, 



043506 



Viel M, 

L51 
Viel M. 

684 
Viel M. 
Viel M 



Haehnelt M. G., Lewis A., 2006, MNRAS, 370, 
Haehnelt M. G., Springel V., 2004, MNRAS, 354, 



Springel V., 2010, JCAP, 6, 15 
Haehneh M. G., Matarrese S. 



Haehnelt M. G.^ 
Lesgourgues J. 
Riotto A., 2005, Phys. Rev. D, 71, 063534 
Viel M., Lesgourgues J., Haehnelt M. G., Matarrese S., 

Riotto A., 2006, Physical Review Letters, 97, 071301 
Wang J., White S. D. M., 2007, MNRAS, 380, 93 
Zavala J., Jing Y. P., Faltenbacher A., Yepes G., Hoffman 
Y., Gottlober S., CatineUa B., 2009, ApJ, 700, 1779 



© 0000 RAS, MNRAS 000, 000-000 



10 




%-10 

1-20 

ft," 

§ -40 



z=5 



(100,512) 
(50,512) 

(25,512) 
(12.5,512) 



-50 
-60 




0.1 1.0 10.0 100.0 1000.0 0.1 1.0 10.0 100.0 1000.0 




10 




z- 


--0 


: 



-10 
-20 








1 ' ': 








-30 








'It ;i 


-40 


- 






'^' ^ 


-50 










-60 








\ ^.\ 



0.1 1.0 10.0 100.0 1000.0 0.1 

k (h/Mpc) 



1.0 10.0 100.0 1000.0 
k (h/Mpc) 



